Ordered bicontinuous double-diamond morphology in subsaturation nuclear matter 
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We propose to identify the new "intermediate" morphology in subsaturation nuclear matter ob- 
served in a recent quantum molecular dynamics simulation with the ordered bicontinuous double- 
diamond structure known in block copolymers. We estimate its energy density by incorporating 
the normalized area-volume relation given in a literature into the nuclear liquid drop model. The 
resulting energy density is higher than the other five known morphologies. 



PACS numbers: 21.65.+f, 26.50.+X, 26.60.-l-c 
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Understanding the form of existence of nuclear matter 
in extreme environments is of importance both from the 
point of view of nuclear many body problem and from the 
context of nuclear astrophysics. Nuclear matter with sub- 
saturation densities (O.lng - ris', rig being the saturation 
density) is believed to exist in the inner crust of neutron 
stars and to appear at stellar collapses. This density 
range corresponds to the transitional region between a 
Coulomb lattice of spherical nuclei and a uniform matter; 
Ravenhall et al. and Hashimoto et al. showed that 
nuclear matter experiences various phases in the course 
of density change. After that, this was elucidated by 
various model calculations 0, 0, H 0, H H . With the 
help of the recent progress of computer power, not only 
the ground state at each density but also the dynamical 
phase transitions between them and the excited states 
were studied by means of first-principle numerical sim- 
ulations with the quantum molecular dynamics (QMD) 
method f iol [111 . Among them, the most basic 

is to understand the form of existence of nuclear matter 
at each density. Watanabe et al. found, in addition 
to the five known morphologies — sphere, cylinder, slab, 
cylindrical hole (tube), and spherical hole (bubble) — , a 
new "intermediate" morphology that is characterized by 
negative Euler characteristic, x < 0- They described this 
as a highly connected spongelike shape, and conjectured 
its relevance to astrophysics. 

However, these five morphologies were known (141] in 
block copolymers of, for example, styrene and isoprene. 
In the nuclear case, the two domains are composed of 
matter and void (or very dilute neutron vapor), whereas 
they are composed of two kinds of polymers such as 
polystyrene (PS) and polyisoprene (PI) in the macro- 
molecule case. Their morphologies are determined by 
a balance between the surface energy and the Coulomb 
energy in the former, whereas by that between the (inter- 
)surface energy and the stretching free energy, like that 
causes rubber elasticity, in the latter. A new morphol- 
ogy was found experimentally in a star block copolymer 
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by Thomas et al. |l5l and in a diblock copolymer by 
Hasegawa et al. UG^ between the PS cylinder and the 
lamella (slab) and between the lamella and the PI cylin- 
der. These experiment determined essentially the correct 
morphology (Fig. 4 in Ref. and Fig. 3 in Ref. T^). 
Soon thereafter the shape of the interface between the 
two microphases was mathematically recognized as the H 
surface — a family of surfaces with constant mean curva- 
ture H 17]. The H surfaces are known to minimize the 
area under the symmetry and volume conservation con- 
dition. The observed morphology is called the ordered 
bicontinuous double-diamond (OBDD) structure accord- 
ing to its symmetry. The OBDD structure consists of two 
interwoven networks of tetrahedral units (four fold junc- 
tions) filled by one material and the remaining matrix 
filled by the other. Calculations of its free energy were 
done by several groups JJt, ,1^ J^J and they concluded 
that the OBDD structure is not the ground state at any 
PS/PI composition. 

Both from the location — adjacent to the slab — of the 
"intermediate" phase of Watanabe et al. and their ob- 
servation that it is highly connected spongelike, it looks 
quite natural to interpret this phase as the OBDD struc- 
ture observed in block copolymers. A direct calculation 
of the energy density of this structure based on some mi- 
croscopic nuclear Hamiltonian is desirable, but it is too 
much complicated unfortunately. Alternatively, here we 
estimate its energy density by combining the familiar liq- 
uid drop model relation and the normalized area-volume 
relation of the OBDD morphology given in the litera- 
ture 

The liquid drop model relations are taken from Raven- 
hall et al. who first predicted the non-spherical mor- 
phologies. Under the Wigner-Seitz cell approximation, 
the total energy density is given by Eg + Ec -\-EB-\-Ee, a 
sum of the surface energy, the Coulomb energy, the bulk 
energy, and the kinetic energy of electron gas. Here the 
Coulomb energy consists of the nuclear electrostatic en- 
ergy and the lattice energy of electron gas and a spatially 
spread nucleus embedded in it. Zero temperature is as- 
sumed. The model is formulated by extending that for 
the spherical case given in Ref. |23| • We consider spheri- 
cal (c? = 3), cylindrical (d = 2), and slab (d = 1) cases as 
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in Ref. . The unit cell for each case is a sphere with the 
radius rc, a cylinder with the radius and the length /, 
and a rectangular parallelepiped with the sides 2rc, a, &, 
respectively. The average density over the cell is n. In 
each cell, nuclear matter is put in the form of a sphere 
with the radius r, a cylinder with the radius r and the 
length and a rectangular parallelepiped with the sides 
2r, a, h, respectively. The density of nuclear matter is 
n' . The other part of the cell is occupied by a very dilute 
neutron vapor. The proton fraction is x. The volume 
fraction is m = n/n' — [r/rcY- The surface tension is a. 
Among the total energy density, and Eq depend on 
the shape and size of the cell and given by 



Es 



uad 



r 

Ec = 27Tn'^x^e\^ufa{u), 
1 r 2 / 1 , 1 2 
d-2\ 2 



fdiu) = 



d + 2 



lira fdiu) 

The variation of Es 
familiar relation 



1 — In u 



Eq with respect to r gives the 



Es = 2Ec. 
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The bulk energy is given by 
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with Eq and K being the binding energy per baryon and 
the incompressibility. The variation of Eg + Ec + En 
with respect to n' gives an equation that determines u [ij . 
That for the cylindrical and spherical hole morphologies 
is similar. Thus, n', r, and Tc for each n and morphology 
are determined. The electron energy that is common to 
all morphologies are irrelevant to energy comparison but 
can be given by 



Ee = —hc(^3TT^ne^ He, 



He = Xn, 

as an ultrarelativistic gas [2^1 . Adopting the parame- 
ter set X = 0.3, Eq = -11.4 MeV, K = 291 MeV, 
ns = 0.147 fm-3, and a = 0.74 MeV/fm^ given in Ref. Q 
and relevant to stellar collapses, the relative energy den- 
sity and the cell size are obtained as in Figs. and El 
Figure ^ indicates the sequential shape change as the 
density change. The cell sizes of the spherical nucleus 
and hole in Fig. [2 are used later. 

In the above model for the five known morphologies, 
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are given independently and 



accordingly the variation of their sum with respect to r 
leads to Eq. QJ. In the case of the OBDD structure, how- 
ever, Eg and Ec can not be represented simply (at least 
to the author's knowledge). Alternatively, we can utilize 
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FIG. 1: (Color online) Total energy density of each morphol- 
ogy relative to that of the slab, as a function of the average 
density. 
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FIG. 2: (Color online) Cell size of each morphology as a func- 
tion of the average density. 



the non-dimensionalized area- volume relation given for a 
cubic cell in Ref. 0. With a lattice parameter X, this 
gives S/X"^ as a function of V/X^ . Consequently, the 
surface energy density is given by 
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as a function of u = V/ X'^ . The relation between u and n 
is taken from the d = 3 (nucleus or hole) case. Assuming 
that Eq. ^ holds also for this morphology, the Coulomb 
energy is automatically determined. This means that 
the total energy can be obtained since E-q and E^ are 
independent of morphology. The area-volume relation 
in Fig. O for the OBDD morphology was adapted from 
Fig. 1(b) in Ref. (171 for the single-diamond structure. 
Since the relation 
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can be derived form the first variational formula of area, 
the nuclear OBDD {V/X^ < 0.5) is a family of surfaces 
with H >0, while the hole OBDD (V/X^ > 0.5) is that 
with H <0 (see Fig. 4 in Ref. [13). Ahhough the H > 
and H < parts are connected smoothly in the case 



of the single-diamond structure, their curvatures are dis- 
continuous in the case of the double-diamond structure. 
This indicates that the lamella structure exists between 
them. 
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FIG. 3: Normalized area-volume relation for the OBDD struc- 
ture. This is adapted from that for the single-diamond struc- 
ture in Ref. p^ . The derivative is proportional to the mean 
curvature. 

Assuming X — l^j rdd — 3) (see Fig. EJ, we esti- 
mate the energy density of the OBDD phase. The result 
is shown in Fig. 0] This figure indicates that this sim- 
ple estimate gives 25 - 30 keV/fm'^ higher energy for the 
OBDD structure than the slab. Qualitatively, this re- 
sult is consistent with that the OBDD phase is not the 
ground state at any composition in block copolymers. In 
the previous works, the "cross" phase in Ref. 4] and the 
"mixed" phase in Ref. |^ might correspond to the OBDD 
phase although our estimate gives higher energy. We did 
not try to change the parameters given in Ref. Q be- 
cause all parameters correlate and they require Skyrme 
model calculation that is beyond the scope of the present 
simple estimate. 

To summarize, we proposed to identify the new "in- 
termediate" morphology in subsaturation nuclear matter 
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FIG. 4: (Color online) Total energy density of nuclear and 
hole OBDD structures relative to that of the slab, as functions 
of the average density. Those for the known morphologies are 
the same as in Fig. 

observed in a recent QMD simulation with the ordered 
bicontinuous double-diamond structure known in block 
copolymers. We estimated its energy density in a hybrid 
manner — incorporating the normalized area-volume re- 
lation given mathematically in a literature into the nu- 
clear liquid drop model. The resulting energy density is 
higher than the other five known morphologies; this is 
qualitatively consistent with the results for block copoly- 
mers. 

Note added in proof 

After submission of the manuscript, the author found 
that another constant mean curvature surface, the bi- 
continuous double-gyroid structure similar to the OBDD 
but consisting of three- fold junctions j21|, is favored in 
block copolymers. In the nuclear case, mathematical sur- 
faces are meaningful as an idealization and therefore it 
would be difficult to distinguish them. 

The author thanks S. Ei for directing his interest to 
block copolymers. 
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